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ABSTRACT 

We  study  the  variational  problem 

2  l«l  2 

j('j')  =  /  |v<i>n  +  /  (4>*‘(v)rdv 

n  o 

i|i  =  0  on  3ft,  =  0,  \|»*(|ft|)  =  1 

where  \ p*  is  the  increasing  rearrangement  of  ty*  An  approximate  problem  is 
introduced  which  involves  a  variational  problem  with  n  free  boundaries 
(n  +  00 ).  Various  estimates  are  established.  In  particular  when  ft  is 
convex  it  is  shown  that  the  solution  to  the  approximate  problem  is 
superharmonic  and  has  bounded  gradient. 
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SIGNIFICANCE  AND  EXPLANATION 


Queer  differential  equations  first  arose  in  the  work  of  Harold  Grad  on 

controlled  thermonuclear  fusion.  They  relate  in  particular  to  models  for  the 

r‘  ' 

slow  adiabatic  evolution  and  resistive  diffusion  of  a  plasma.  They  are  queer 
in  that  they-  share  aspects  of  partial,  ordinary  and  functional  differential 
equations.  In -the -present  work  the  authors  give  a  new  way  of  thinking  of 
these  equations  by  relating  them  to  free  boundary  problems.  This  is  the  first 
of  a  series  of  papers  intending  to  demonstrate  that  solutions  of  such  a  queer 
differential  equation  can  be  thought  of  as  limits  of  solutions  of  free 
boundary  problems  with  n-free  boundaries.  A  long  term  hope  is  that  this  work 
will  complement  and  further  refine  existing  numerical  schemes  for  finding  such 
solutions.  (  - 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC ,  and  not  with  the  authors  of  this  report. 


A  NEW  APPROACH  TO  "QUEER"  DIFFERENTIAL  EQUATIONS 


Peter  Laurence  and  E. 


W.  Stredulinsky 


Introduction . 

In  this  work  we  study  the  so  called  queer  differential  equations.  We  introduce  a  new 
approach  which  gives  insight  into  the  types  of  estimates  one  can  expect  and  which  provides 
a  setting  in  which  it  is  hoped  they  can  be  rigorously  established.  In  this  approach  we 
introduce  an  approximate  problem  involving  many  free  boundaries.  Exploiting  the 
associated  free  boundary  conditions  and  properties  of  harmonic  functions  we  show  that 
solutions  are  superharmonic  and  satisfy  gradient  bounds.  It  appears  likely  that  these 
estimates  will  carry  over  to  the  original  problem. 

Queer  differential  equations  (Q.D.E.s)  arise  in  nuclear  fusion  research.  They  were 
introduced  by  Harold  Grad  as  a  way  of  more  accurately  modeling  adiabatic  compression  [8] 
and  resistive  diffusion  [5]  of  a  plasma. 

Historically  "plasma  equations"  have  been  studied  by  many  authors.  These  studies 
have  essentially  involved  nonlinear  elliptic  equations  and  the  associated  free  boundary 
problems.  Other  models  which  lead  to  Q.D.E.'s  have  been  considered  by  J.  Mossino, 

R.  Temam  [10.11),  G.  Vigfusson  [15,16].  The  origin  of  all  the  plasma  equations  is  the 
Grad-Shaf ranov  equation  which  in  2-dimensions  (the  case  of  physical  interest)  is 
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Here  f(V)  ere  profile  functions  related  respectively  to  the  pressure  and  the  total 


poloidal  current  [6].  Traditionally  p,  f  have  been  prescribed. 

The  work  of  Harold  Grad  demonstrates  that  the  quantities  which  are  prescribed  by  the 
dynamics  of  the  plasma  are  not  p ('If),  f(^)  but  rather  u(4>)  end  v(\Ji)  given 
respectively  by  (see  [5]) 

(2)  U<V)  -  - E— 

<**■ )T 


where  y  *  2  in  the  two  dimensional  case.  Here  iji*  is  the  increasing  rearrangement  of  a 
function.  To  recall  its  definition  we  first  introduce 

(3)  V^(t)  -  |(x  :  *(x)  <  t} | 

(where  |e|  is  the  Lebesgue  measure  of  E).  Then  i(i*(v)  is  the  inverse  function  of 
V^.  More  exactly 

(4)  ♦*(v)  «  inf{s  s  V^(s)  >  v}  . 

Eliminating  p,  f  from  (1)  in  favor  of  u  and  v  leads  to 

(5)  A*  -  -|l*  <♦)  <+**  )Y  - 

-  ^  Cv2(H»)>*  (♦**  )2  -  v2(*)i|»*"  . 

Here  4>*'  and  ♦*"  denote  the  compositions 

♦*'  <V+(#(x)  > )  and  **"(VH|(*(x))) 

respectively.  From  (5)  it  is  clear  that  in  a  typical  case  second  derivatives  of  i|>* 
appear.  For  this  reason  we  propose  to  study  the  model  equation 


(6) 

A*(x) 

“  -**"<V<,(*(x))) 

obtained  in  2-dimensions  (where 

Y  *  2) 

from 

the  choice 

-  V2  , 

V(ljr)  -  0  . 

An  equation  of  this  form  was  first  considered  by  R.  Temam  in  (13)  where  a  variational 
formulation  was  given  and  the  existence  of  a  minimizer  was  established.  There  the  author 
considered  a  domain  £J  C  R°  and  imposed  the  condition  ♦  ”  0  on  3fl.  However  physically 

(6]  what  one  would  like  to  also  prescribe  is  the  difference  in  tf  (the  flux)  between  its 


absolute  maximum  and  minimum  (here  the  minimum  is  zero).  This  is  equivalent  to  specifying 
^*(0)  and  i|»*<  We  will  take 

(7)  **(0)  *>  0,  **(  |ft|)  -  1  . 

This  is  no  loss  of  generality  as  it  is  easily  checked  that  (6)  is  invariant  under  the 
transformation 

♦  *  l  (i  +  c 

for  constants  A,  c. 

In  section  1  of  this  paper  we  will  begin  as  in  [13J  by  formulating  a  variational 
principle  in  Q  (an  open  set  in  if1)  and  establishing  the  existence  of  a  minimizer  in  a 
function  class  which  incorporates  (7).  we  then  derive  certain  estimates  for  the 
minimizer.  Obstacles  towards  obtaining  further  estimates  are  the  apriori  unknown 
structure  of  the  level  seta  and  of  the  set  of  critical  points.  Thus,  as  an  attempt  to 
isolate  regularity  questions  from  these  we  introduce  in  section  2  a  variational  problem 
within  a  class  of  functions  with  convex  level  sets.  In  section  3  we  approximate  this 
problem  by  an  n-ahell  free  boundary  problem  which  makes  clearer  the  separate  influences  of 
the  Laplacian  term  and  the  Q.D.E.  term.  Similar  problems  with  one  free  boundary  were 
first  studied  by  Acker  [1]  and  Caffarelli  and  Spruck  [33.  In  a  future  paper  we  intend  to 
show  that  a  solution  of  (6)  is  obtained  as  a  limit  of  solutions  to  the  approximate  problem 
as  n  ♦  •• 

In  the  present  paper  in  section  4  we  show  that  the  solution  to  the  approximate 
problem  is  superharmonic  and  satisfies  a  gradient  bound.  Superharmonicity  is  preserved 
under  taking  a  weak  limit  of  such  solutions.  There  is  strong  evidence  that  such  a  limit 
is  a  solution  of  the  Q.D.E.  (6)  and  that  the  gradient  bounds  are  preserved  although  all 
details  have  not  been  worked  out  as  of  this  writing.  Recent  work  of  Caffarelli  and 
Friedman  [2]  may  make  possible  the  extension  of  this  approach  to  the  nonconvex  case. 

We  would  like  to  thank  H.  Grad  for  suggesting  this  problem  and  for  his  warm 
encouragement  and  helpful  conversations  throughout  this  work.  In  addition  we  would  like 
to  thank  L.  Caffarelli,  A.  Friedman,  J.  Spruck  and  W.  Ziemer  for  useful  conversations. 
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1 .  A  Variational  Problem  and  Estimates  for  Its  Mlnlmlzers 

In  this  section  we  will  introduce  a  variational  formulation  for  the  problem  (6). 

This  is  similar  to  the  variational  problem  considered  in  [13],  however  in  [13]  boundary 
conditions  were  not  imposed  on  i]i*.  Vie  begin  with  a  brief  discussion  of  the  formal 
relationship  between  the  Q.D.E.  (6)  and  the  variational  problem  (8,9).  We  then  establish 
the  existence  of  minimizers  and  go  on  to  prove  various  estimates  for  stationary  solutions 
of  the  variational  problem.  For  instance  it  is  shown  that  1*  Lipshits. 

Let  Q  be  an  open  bounded  set  in  VP  and  define 

,  2  J°l  2 

(8)  J(i[i)  *  /  !Vf|  dx  +  /  (♦*')  dv  . 

ft  0 

We  seek  to  minimize  J(^i)  over  in 

(9)  w  -  {+  e  w1,2(fl)  :  e  w1,2<o,|ft|),  **<o)  -  o,  v*|n|)  -  i}  . 

By  w1,2(E)  we  mean  the  Sobolev  space  of  square  integral  functions  on  E  with  square 
integral  first  derivatives  and  by  W1,2(E)  we  mean  the  closure  of  c“(E)  in  W1,2(E>. 
Clearly  we  have  imposed  the  boundary  condition 

i|i  »  0  on  3ft. 

In  addition  the  constraints  <M0)  »  0,  i)i*(|a|)  -  1  are  equivalent  to  specifying  the  inf 
and  sup  of  i|i  to  be  zero  and  one  respectively.  The  condition  ♦  (0)  “  0  is  actually 
superfluous,  see  Remark  1.1. 

We  will  now  give  a  brief  derivation  of  the  Euler  equation  for  (8,9).  The 
calculations  are  somewhat  formal  but  illustrative. 

2  — 

Suppose  t|>  is  a  minimizer  or  simply  a  stationary  point  and  i|>  e  C  (ft).  Let 

pfc  -  {x  e  ft:  <»(x)  »  t}. 

Since  e  C2(ft)  we  have  that  Q\(p0  u  Pi>  i*  °P®n  an<5  consequently 

sup  $  <  1,  inf  i|i  >  0 
*  K  -  „ 

for  any  compact  set  K  C  fl\(pg  U  P|).  Given  1  8  Cg(B\(pp  u  Pi))  we  then  have 
^•^t  +  enew  for  small  e  and  we  can  carry  out  the  following  variation. 


(10) 


0  -fe  J(V)e-0  *  2  /  .  vn  dx 


+  2  /  1(1* ’  <  n  >'dv 

0 


where 


with 


*(v)  -  {x  :  V^(<p(x))  -  v> 

and  Hn_1  being  n  -  1  dimensional  surface  measure  (Hausdorff  measure).  The  second  term 

on  the  right  hand  side  of  (10)  Is  obtained  by  using  the  following  formula  along  with  the 

assumption  that  differentiation  with  respect  to  e,  v  can  be  interchanged, 

+  en)*(v)  -  <  n  >  (v)  . 

'e-0  _ 

This  formula  can  be  rigorously  derived  under  the  assumption  that  |Vi|/|  ft  0  in  0,  [11]. 

An  integration  by  parts  now  yields. 


/  nAi>  dx 

a 


lot 
♦  J 


l|l*"  <  n  >  dV  m  0 


which  can  be  expressed  after  use  of  the  coarea  formula  (12)  as 

/  ni*  dx  +  /  ♦*"(V  <*(x)))n(x>  dx  -  o  . 

n  q  v 


Therefore 

A<|>  -  in  £i  \(Pq  UP,). 

It  is  possible  to  perform  variations  involving  n's  whose  support  does  not  He  in 
P0  U  p,  (see  remark  1.4).  These  however  bring  into  play  certain  boundary  terms  due  to 
the  integration  by  parts  in  the  second  term  of  (10).  The  interpretation  of  such  terms 
requires  an  understanding  of  the  set  of  critical  points  of  the  minimizer  y  and  its  level 
sets . 
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He  now  present  an  existence  result  and  various  estimates  for  stationary  points  of  the 
variational  problem.  For  reasons  of  clarity  we  defer  the  proofs  until  later. 

Theorem  1.1.  The  infimum  of  J(<0  in  W  is  attained. 

Remark  1.1.  There  is  a  weak  maximum  principle  for  the  variational  problem  in  the 
following  sense.  It  is  not  necessary  to  impose  the  condition  ( 0 )  •  0.  If  this 
assumption  is  dropped  from  the  definition  of  H,  theorem  t.1  remains  true  and  it  can  be 
demonstrated  that  ip*(0)  «  0  still  holds  for  any  minimizer  See  the  end  of  the  proof 

of  theorem  1.1  for  details. 

Remark  1.2.  It  is  not  known  whether  the  minimizer  is  unique  in  general.  It  is  unique 
however  in  the  radial  case,  see  remark  I.S. 

He  would  now  like  to  introduce  a  regularity  result  for  To  do  this  we  introduce 

a  function 

W(t)  -  /  |V<r|  2dx  . 

{^»<t} 

Clearly  H  is  an  increasing  function.  Let  dH  be  the  induced  Lebesque-Stielt jes 

measure.  Note  that  if  we  interpret  i|>  properly  (see  (12))  it  follows  from  the  coarea 

formula  that  dH  is  absolutely  continuous  and 

^f(t)  •  /  | Vi|» | dHn” 1  a.e.  . 

dt  {*-t> 

Theorem  1.2.  If  lji  is  a  stationary  point  for  the  variational  problem  (8,9)  then 

H  +  ♦•'(V^)  -  J(*)  a.e. 

and  therefore  i)»*  is  Lipshitz. 

Remark  1.3.  The  methods  of  theorem  1.2  apply  to  the  problem  of  minimizing  J(4>)  over 

♦  e  H  having  a  given  collection  of  level  sets  “  t}  -  H(t).  This  is  reminiscent  of  a 

numerical  scheme  introduced  by  Harold  Grad  [7],  see  remark  1.6. 

In  order  to  present  several  further  estimates  we  need  a  general  version  of  the  coarea 

1  2 

formula.  In  this  respect  we  discuss  averages  of  Sobolev  functions.  For  fen'  (SI)  let 


f 

r 


(x) 


/ 

Br<x) 


f(y)dy 


where  Br(x)  =  {y  :  | x  —  y |  <  r}.  Also  let  X(x),  u(x)  be  the  approximate  lower  and 

upper  limits  respectively  of  f  at  x  (see  [4],  2.9.12)  and  X(x)  -  liminf  fr(x), 

1  2  r>0 

U(x)  «  limsup  f  _  ( x )  •  It  is  easily  shown  that  f  ♦  f  in  M  '  (fl).  From  this  it  is  seen 
r*0  r  loc 

that  f  ♦  f  pointwise  quasi-everywhere,  that  is  everywhere  except  possibly  on  a  set  of 

zero  Newtonian  capacity.  Let  f(x)  **  Uit  f  (x)  where  it  exists.  Clearly 
12  ^0 

f  “  f  Xn  w  '  (A).  Throughout  the  rest  of  this  section  the  particular  representative 
1  2 

(from  aw'*  equivalence  class)  which  we  use  will  be  the  above  limit  of  averages  (which 

is  defined  quasi-everywhere).  Using  the  fact  that  f  ♦  f  in  w’,2(ft)  it  can  be  shown 

r  loc 

in  addition  that  X(x)  =  X(x)  »  u(x)  «  y(x)  quasi-everywhere.  Using  this,  part  14  of 
theorem  4.5.9  [7]  and  the  fact  that  the  Newtonian  capacity  of  a  set  being  zero  implies 
that  its  Hn_1  measure  is  zero  we  get  the  following  result, 


(12)  /  g|Vf|  -  j  j  g  dHn-1dt 

n  —  {f-t> 

for  measurable  functions  f,g  such  that  g|Vf|  is  integrable  and  few1'2  (with  the 
particular  representative  of  f  being  chosen  aB  above). 

Using  (12)  we  can  now  establish  the  following  lemma. 

Lemma  1.1.  Let  dvu(t)  be  the  Lebesgue  Stieltjes  measure  induced  by  Vu(t)  for 
u  e  w1,2((l)  then 

i)  dVu(t)  -  (u*'(Vu(t)))_1dt  +  l  |{u  -  t  )td6  (t) 

i  i 

where  d6t  (t)  is  a  unit  mass  concentrated  at  tA  and  { t is  the  collection  of 
points  t  for  which  |{u  ■  t}|  >  0. 

ii)  dV  <t)  -  /  |Vur1dHn-1dt  +  ds(t) 

(u-t) 


where  S(t)  »  | {u  <  t)  n  { | Vu |  -  0 } | 


Proof  of  Remark  1 . 1  (weak  maximum  principle).  Remove  the  condition  i(i*(0)  «  0  from  the 
definition  of  W.  The  proof  of  theorem  1.1  still  holds.  Assume  4>  is  a  minimizer.  let 

-  Kax{i|>,0}  . 

°  1  2 

It  is  easily  checked  that  (Y+)*  *  +  e  Also  i|i+  6  M  '  (ft)  with 

7*+  “  X^o}7*  and  6  W1,2(0,|fl|>  with  •  Here 

XE  “  1  on  E  and  xR  "  0  on  the  complement  of  E.  Clearly  <  J ( vp )  unless 

+  >  0  a.e. 

Remark  1.4.  This  weak  maximum  principle  in  conjunction  with  a  rescaling  leads  to  the  fact 

that  W  may  be  replaced  by 

_  O  .  «  i  2 

(13)  w  -  e  w  '  (Si)  i  **  e  w  '  (o, |o| >,  **i ja| >  > 

without  altering  the  infimum  of  J. 

Remark  1.5.  In  the  case  when  0  is  a  ball  the  variational  problem  has  only  one  radial 
stationary  point.  This  is  in  fact  the  unique  minimizer.  To  see  this  let 
0  «  (x  i  | x |  <  R}  and  assume  t|i(x)  “  f(|x|)  is  a  stationary  point.  The  problem  is 
easily  reduced  to  considering  an  appropriate  O.D.E.,  for  which  there  is  a  unique  solution 
satisfying 

f’(r)  -  -J<*)(nw  rn_1  +  (nu>  r""1)"1)"1 
n  n 

where  is  the  measure  of  the  unit  ball  in  S^.  By  using  Bymmetrization  (which  leaves 

t|i*  fixed)  we  see  that  any  minimizer  must  be  radial  and  therefore  unique. 

Proof  of  Theorem  1.2.  Let  “  f£o<(i  where  f£(x)  «  x  +  eh(x),  h  e  c”(0,1)  and  e  is 
small  enough  that  eh’  <1  on  [0,1].  We  have 


?<t>e  “  f^(¥>V*  . 

00 

Also  f  is  strictly  Increasing  with  C  inverse  so 


Consequently 


and 


Wf(t,)  ■  vtl  • 

(  f  Olp  >  *  (  V  )  «  f  (l|»*(v)  ) 


( f oV) * '  ( v)  -  f ' (¥*(v>)**’<v)  . 
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From  this  it  is  seen  that 


2  2  11  2  2 
J(#_>  -  /  (f'<*))  |V*|  d x  +  /  (C(**(v))  (♦••(v))  dv  . 

e  a  C  0  e 


It  is  easily  seen  that  the  following  formula  holds  for  continuous  g. 


/  g(t)dWft)  -  /  g(4t)|V*[2dx  . 

o  n 


Therefore 


2 

/  h* (t)dw(t)  -  /  h'(<>)|V*|  dx  . 

0  a 

Furthermore  if  r(v)  is  increasing  and  absolutely  continuous  and  g  is  Borel  measurable 
then  (see  [12]) 


b  r(b) 

/  g(r(v) )r* (v)dt  »  /  g(t)dt  . 
a  r(a) 


To  apply  this  formula  we  introduce  the  following  device*  Let 


where  V^t*),  V^ft") 
We  then  have 


%<t) 


V*  >  *  V*~> 


represent  the  limit  of  from  the  right  and  left  respectively. 


**’ (v)  -  ♦*' CV^(**( v) ) )  a.e. 

since  off  of  a  countable  set  either  both  sides  are  zero  or 

v  -  V^(<i»(v))  . 

Consequently  letting  r(v)  -  <i*(v),  g(t)  -  h' (t)^*' (v^(t)  )  in  this  formula  we  get 
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Therefore 


1 

-  I 

0 


i 

0  -  /  h' (t) (dW(t )  +  ^*'(V^(t))dt) 
0  v 


and 

(14)  dW(t)  +  <f*' (V^(t))dt  -  c  a.e.  . 

Clearly  if*  Is  Lipshits  since  dw(t)  >  0.  Integrating  (14)  over  (0,1)  and  reversing  the 
above  arguments  leads  to 

c  -  J(if). 

Remark  1.6.  The  methods  used  above  are  closely  related  to  a  more  geometrical  way  of 
looking  at  the  variational  problem.  This  approach  which  we  outline  below  is  made  use  of 
in  numerical  schemes  devised  by  Harold  Grad.  The  relationship  between  it  and  the 
variational  problem  (6,9)  will  remain  formal  until  it  can  be  proved  that  there  are  no 
Intervals  on  which  </*'  *  0  when  f  is  a  minimi zer  to  (8,9). 

Let 

V(x)  -  ( if  ( x ) )  . 

Clearly  V  is  determined  solely  by  the  level  sets  of  if  and  thus  is  a  purely  geometrical 
quantity.  Let 

K(v)  -  _  J  | Vv | dHn” 1  . 

{V(x)-v} 


It  is  possible  to  show  that  <7(ii)  is  minimized  by  if  satisfying 


lf.l(v)  - 


1  +  K(v) 


with  J(< f)  being 
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1  131 


l  J 
0 


1  +  K(v) 


U  Y  J 


The  problem  of  minimizing  (8)  can  thus  be  viewed  as  finding  a  collection  of  level  sets  for 
which  (IS)  is  minimized. 

Proof  of  Lemma  1.1.  By  [14]  u*  e  w1,2(0,a)  for  any  a  <  |&|.  We  establish  the  result 
on  (0,a)  from  which  the  general  result  follows. 

Remove  the  "flat  places"  in  the  graph  of  u.  ltoat  is  let  {(an»bn)}  the 
collection  of  intervals  on  which  u*  is  constant  and  define 

u*(v)  ”  u*(v) 

for  v  m  v  -  \  (bn  -  an)  and  v  8  U(an,bn).  u*  is  absolutely  continuous  and  strictly 


b  <v 


increasing  so  ?t  has  absolutely  continuous  inverse 


g  with 


g’(t)  -  (u*Mg(t)))-1  a.e.  . 


Also 

u**(v)  -  u** (v) 

for  almost  all  v  e  (0,a) \ (^(an,bn) ) ,  in  fact  precisely  for  those  v  which  are  Lebesgue 
points  for  the  set  (0 ,a)\ (^(an ,bn) ) .  It  is  easy  to  reconstruct  Vu(t)  as  the  sum  of 
g  and  an  increasing  step  function  with  jumps  corresponding  to  the  "flat  places"  on  the 
graph  of  u  and  so  i)  follows. 

We  can  conclude  ii)  from  an  application  of  the  coarea  formula.  To  see  this  let 

9  “  XE(tJ7ul  >  Of  “  < 

so  that 

|*(t)|  «  /  g|Vu|dx 

a 


■  /'  I  dH"-' 

minu  (u“s) 


TvuT 


ds 


But 


0  “  l  *{|Vu|-0}!Vuldx  "  / 


lit) 


■0} 


dHn_1dt 
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7V*'V"5r’.%T%T%-' 


80 

(16) 

therefore 


Hn  ^ ( {u  •  t,  |Vu|  *  0})  *  0  a.e. 


vu(t)  -  /  /  |7ui-1dHn-1  ds  +  |(|Vu|  -  0,  u  <  t}( 

minu  {u”s} 

and  ii)  follows. 

Since  by  li)  ds(t)  contains  the  singular  part  of  dvu(t),  that  is  (by  i) ) 

I  |{u  *  t.} |dfi  (t)  , 

1  i 

it  is  clear  that  iii)  must  hold. 

Proof  of  Proposition  1.1.  Let  s(t)  ”  (4>“t) .  For  almost  all  t  we  have 

/  |7*|  <  |V*|  -  0}>  -  0 

s(  t ) 


(see  (16))  and  lemma  1  iii)  holding  with  u  *  i|i.  For  such  a  t  choose  e  such  that 

I  |V*|dHn-1  -  EHn_1(s(t)) 

8(t) 

E  ■  {xes(t)  «  |V(|»|  <  2e}  . 


and  let 


From  this  we  see  that 


but 


2eHn"’(s(t)\E)  <  eHn-,(s(t) ) 


Rn-,(s(t))  <  2Hn-1(E) 


**'(v  <t))  <  (  /  Ivvl-V  <  2e(Hn'1(E))-1 
V  8(t) 

<  4E(Hn_1(s(t)))'1 


4 (Hn_1  (  s(  t ) ) ) ”*  |  |Vil>|dH  • 

s(t) 


Combining  this  with  theorem  1.2  gives  i). 
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We  now  wish  to  use  the  following  "i soper imetric  inequality*, 

n- 1 

(17)  (|(1|  -  V  (t))  n  <  iT’u  “1/nHn~1((t  -  t})  a.e. 

y  n 

where  u>  is  the  volume  of  the  unit  ball  in  R*1.  For  smooth  a  it  is  a  direct 
n  T 

consequence  of  the  standard  isoperimetric  inequality  and  the  fact  that  almost  all  level 

r2 

sets  of  41  are  smooth  (Sard’s  theorem,  implicit  function  theorem).  For  ♦  e  W  '  (8)  one 
can  combine  the  Sobolev  inequality 

n- 1 

{  /  un/<n_1)}  n  <  n'V1/n  /  | Vu | 
n  n  0 

with  the  coarea  formula  (12)  to  prove  (17). 

By  (16)  we  have 

Hn  1  ( {41  «  t,  |y*i|  -  0)  -  0  a.e. 
so 

(Hn-1({*  -  t}))2  <  /  IVtf’dH"-1  j  |V*|dHn_1  a.e. 

{♦-t}  t*«t} 

»  /  iv»r1dHn-1(j(*)  -  ***(v .  (t) > » 

{*-t}  * 

by  theorem  1.2.  Now  combining  this  with  leimna  1.1  iii)  and  (17)  we  get 

2(n-1)  -1 

**<w <  j<*’ii  +  ('°i  -  vtn  "  i 

almost  everywhere.  We  can  easily  recover  ii)  since  is  sero  on  (a^,b^),  the 

intervals  where  <t>*  is  constant,  and  if  F  has  measure  sero  then 

|V<I((0,1)\F)|  -  |(0,  |a|)  |  -  l  (bA  -  8i) 

(recall  the  proof  of  lemma  1.1  i)).  From  remark  1.5  it  can  be  seen  that  ii)  is  sharp  if 


0  is  a  ball. 


2.  The  Variational  Problem  for  Functions  with  Convex  Level  Sets 

For  a  convex  domain  fl  we  expect  minimizers  to  the  variational  problem  (8,9)  to  have 

nested  convex  level  sets.  This  motivates  us  at  the  present  stage  to  introduce  such  a 

condition  as  a  constraint  on  the  variational  problem.  The  hope  is  that  solutions  of  the 

constrained  problem  will  correspond  to  solutions  of  the  original  problem.  This  is 

justified  to  some  extent  by  the  results  to  be  presented  in  [9] .  For  simplicity  we 

2 

restrict  our  attention  to  ft  C  *  •  This  is  the  case  of  physical  interest. 

He  will  first  present  the  technical  setting  for  the  constrained  problem  and  then 
prove  existence  and  continuity  of  solutions. 

Let  us  begin  by  giving  a  weak  definition  for  the  class  of  functions  with  convex  level 
sets.  He  have  chosen  this  particular  definition  since  it  accomodates  the  "almost 
everywhere  defined”  nature  of  functions  in  H1,P(Q)  and  more  importantly  allows  us  to 
avoid  the  technicalities  associated  with  "flat  places”  in  the  graph  of  a  function,  that  is 
places  where  level  sets  have  positive  measure. 

Por  an  open  convex  set  0  let 

(18)  Cp  •  (u  e  W1,p(fi)  i  There  is  a  dense  set  {tn>,  tn  e  (inf  u,  sup  u)  and  there 

are  convex  sets  Cn  such  that  u  >  tn  a.e.  in  Cn  and 

u  <  t  a.e.  in  fl\C  }  . 
n  n 

In  this  context  it  is  easy  to  see  what  form  "flat  places"  take.  To  see  this  consider  a 

function  u  8  C_.  Let  {t  ),  {C _}  be  as  in  (18).  Given  t  €  (inf  u,  sup  u)  choose 
p  n  n 


subsequences  {t,  },  {t,  }  of 

» t  n  2  $  n 

U  } 

n 

such  that  t,  +  t 

i  ,n 

and 

fc> 

3 

+  t 

and  let  the 

associated  convex  sets  be  (n. 

c2,n* 

It  is  easy  to  check 

that 

{u  - 

t} 

differs  from 

(19) 

Ht)  - 

<Q  C,  >\(U  c  ) 

n  i  ,n  n  2  #n 

by  a  set  of  measure  zero.  If  | {u"t} {  >  0  then  T(t)  is  a  region  bounded  by  the  convex 


C.  ),  3(U  C  )  . 

n  i  (ii  n  z  f  n 

If  |{u”t}l  »  0  then  f(t)  is  a  convex  curve.  In  the  case  where  t  «  inf  u  or 


t  *  sup  u  one  can  use  c2,n’  c1,n  respectively  to  show  that  {u»t}  differs  from 
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%  . 


5 


r.- 

!■  . 


inf { j(t )  S  *  e  w  nc2>  . 


Q\<n  c2,n^  and  n  c1,n  respectively  by  a  set  of  measure  zero.  In  the  first  case  a 

"flat  place"  is  bounded  by  3ft  and  3  (^  C,  )  and  in  the  second  by  3  (*"»  C,  _ ) • 

n  4,n  n  •#>» 

we  can  now  state  the  constrained  problem  as  that  of  minimising  J(f)  over  Oj, 

that  is 
(20) 

Theorem  2.1.  The  infimum  in  (20)  is  attained  in  W  nC2« 

Remark .  Theorem  2.1  is  an  immediate  consequence  of  the  proof  of  theorem  1.1  and  the 
following  result  on  the  weak  closedness  of  C2* 

Proposition  2.1.  Given  e  C  such  that  4n  ♦  4  weakly  in  W1,*>(Q),  p  >  1  and 

t  e  (inf  4»,  sup  40  then  there  exists  a  convex  set  C  such  that  <r  >  t  a.e.  in  C  and 

V  <  t  a.e.  in  ft\C. 

Corollary  2.1.  Cp  is  weakly  closed  in  W^'p(ft)  for  p  >  1. 

In  proving  proposition  2.1  we  need  the  following  lemma.  It  is  actually  true  that 
u*  e  *J£(0.|0|)  (see  [13])  but  only  the  Maker  result  is  needed. 


Lemma 


2.1 .  If  u  e  W1,p(n),  p  >  1  then  u*  is  continuous  on  (0,|ft|). 


In  addition  we  will  prove  the  following  regularity  result. 

Theorem  2.2.  If  u  e  C  ^  n  l" ( ft )  then  u  is  continuous. 

Remark.  Since  Q  C  R2,  w1,p(0)  consists  of  continuous  functions  if  p  >  2. 

Proof  of  Lesma  2.1.  Assume  u*  is  discontinuous  at  v  e  ( 0 , | ft | ) .  There  then  exists 

t2  such  that  inf  u  <  t1  <  t2  <  sup  u  and  Itt^  *  u  <  tj)  |  ”  °* 

1  t  >  t. 


f(t)  - 


*2  -  S 


t,  «  t  <  t2 


t  <  t. 


so  f(u)  e  W1,p(ft)  and  Vf(u)  “  f  •  (u)Vu  a.e.  .  Since  |  C t ^  <  u  <  tj}  |  “  0  we  have 
Vf(u)  ■  0  a.e.  and  also  f(u)  «  f(u)  “  X{u>t  j  a.e.  (xp  equals  1  on  F  and  se: 
0\F).  Since  0  <  v  <  |ft|  we  have  |{u  >  t^}  |  >  0  and  |  {u  <  tj)  |  >  0.  Let 
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E  «  {u  <  tj}.  Using  tha  Sobolev  inequality  (for  |e|  >  0) 

J  U  -  Trr  i  ♦ I p  <  c(o, |e| )  /  |v*|p 
a  |E|  e  n 

with  $  «  f(u)  we  conclude  that 

|{u  >  t2}|  -  0 

which  contradicts  ! (u  >  t  }|  >  0  and  we  are  done. 

2 

Proof  of  Proposition  2.1.  Assume  weakly  in  w1,p(fl)  and  Vn  ®  P  >  1. 

2 

Choose  a  subsequence  (still  represented  by  i|i  )  for  which  tii  ♦  i(i  both  in  L  (Si)  and 

n  n 

pointwise  a.e.  and  iji*  ♦  t|>*  pointwise  a.e..  Given  t  e  (inf  u,  sup  u)  we  can  conclude 
n 

by  lemma  2. 1  and  the  intermediate  value  theorem  that  there  exists  a  v  such  that 

♦*(v)  «  t.  Since  ♦*,  are  increasing  there  is  a  v  <  |ft|  such  that 

ip*(v)  ♦  <Ji*(v)  >  t.  Choose  t_,  C_  from  the  definition  of  *  8  C  such  that  t  ♦  t 

n  **  n  n  p  n 

and  choose  v„  such  that  v „  is  the  largest  v  satisfying  e*(v)  -  t  .  From  the 
»  n  n  n 

monotonicity  of  we  have  vn  <  ▼  for  large  n  so 

lcn I  >  |«|  -  vn  >  |fl|  -  v  >  0  . 

Since  Cn  C  fl.  SI  is  bounded  and  C„  is  convex  there  exist  balls  Bn  C_  Cn  with 

| B^ |  >  X  >  0,  X  independently  of  n.  By  choosing  a  subsequence  and  redefining  the 

balls  Bn  we  can  assume  in  addition  that  they  have  a  common  center.  Parameterizing  3Cn 

using  an  angle  0  we  can  apply  Arzela-ABColi  to  conclude  that  a  subsequence  of  {3  C^} 
converges  uniformly  to  3c  for  some  convex  set  C.  The  desired  conclusion  now  follows 
easily  from  the  properties  of  C^  tn  (see  (18))  and  the  fact  that  ^»n  ♦  i|i  pointwise 
a.e..  The  corollary  follows  immediately. 

Proof  of  Theorem  2.2.  Assume  that  t|i  e  Cj  and  that  we  have  real  numbers  t1f  t2  and 
convex  sets  C1#  C2  such  that  inf  <  t^  <  tj  <  sup  y,  <|>  >  t^  a.e.  on  C^  and 
i(i  <  t^  a.e.  in  £i\C^.  The  following  calculation  will  show  that  3Cj,  3C^  do  not 
intersect. 

Let  d  »  dist  {3c1#  3C2).  Choose  xQ  e  3C^  such  that  d  »  diet  {xQ,3c2}.  Let 
C(r)  be  the  circle  concentric  to  x0  of  radius  r.  Clearly  C(r)  intersects 
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{ i  >  1 2)  and  {41  «  t2)  on  a  a«t  of  poaitiva  ona-dinanaional  measure  for  almost  avary 

r  >  d.  Also  for  almost  all  r,(<  restricted  to  C(r)  Is  an  absolutely  continuous 

function  with  (1-dimensional)  derivative  given  by  T  •  Viji  where  T  is  the  unit  tangent 

to  C(r).  Because  of  this  we  have  for  almost  all  r  >  d  that 

t2  -  t,  -  /  T  •  ViJidH1 

A(r ) 

for  some  arc  A(r)  C  c(r)  n  Si.  Therefore 

(t3  ‘  V*  3  1 

— - - -  *  <=  S  |V*|  2<SH1 

C(r  )no 

and  integrating  with  respect  to  r  from  d  to  1  we  have 

(21)  <t2  -  t,)2  log  d-1  <  c  /  |V*|2  . 

a 

Equivalently  we  have 

(22)  exp(-(t2  -  t,)-2  /  |V*|2)  <  d  . 

Q 

He  will  now  redefine  on  a  set  of  measure  sero  so  that  its  level  sets  {t|>  -  t} 
coincide  with  T(t)  (see(19)).  (21)  then  becomes  an  estimate  of  the  audulus  of  continuity 
of  and  we  are  done  . 

bet  tn,  Cn  he  as  in  the  definition  of  ♦  e  C2«  Let 

r.  -  {x  e  c  «  t  <  t),  r_  -  {x  e  a\c  :  *>t)  . 

n  n  n  n 

Clearly  if  G  -  U  (E„  U  En)  then  |g|  -  0.  Let 
n 

■up  {t  s  x  e  C  }  for  x  e  G 
n  n 

<M  x )  «  ■ 

♦(x)  otherwise  . 

It  is  easily  checked  that  redefining  a.e.  in  this  way  forces  {<i  *  t}  •  T(t)  as 
required. 
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3.  An  Approximating  Free  Boundary  Problem 


We  now  introduce  an  approximate  problem  which  we  will  demonstrate  to  be  a  free 
boundary  problem.  This  is  done  for  several  reasons.  One,  it  regularizes  the  problem  thus 
allowing  estimates  to  be  made.  Two,  it  is  a  first  step  in  achieving  one  of  the  goals  of 
section  2,  that  is  to  produce  a  solution  to  the  original  variational  problem  *  .ich  has 
convex  level  sets.  A  third  advantage  to  this  approach  is  that  it  allows  us  to  seperate 
the  influence  of  the  two  terms  in  J(i(i). 

In  the  next  section  we  will  establish  free  boundary  conditions  and  present 
applications. 

The  approximation  of  J(</)  is  achieved  by  replacing  V*'  by  a  difference 
quotient.  This  yields  the  following  functional.  Let 


(23) 


J_(f)  -  /  l**|2  +  I  fr4 

n  n  i-1  'V 


where  Aa  -  (t1_1  <  Y  <  tj),  tt  *  i/n,  l  «  1,...,7T.  Clearly  f  |  -  “  Vi),(ti-i>' 

We  seek  to  minimize  over 

n 

(24)  r0  -  {*  e  w’''(fl)  n  c2  .  *  (0)  •  o,  ♦  (|n|)  -  i}  . 

It  turns  out  that  J_  does  not  achieve  its  minimum  on  FQ.  What  one  would  expect  to  be 
n 

the  minlmizer  will  be  seen  to  satisfy  <|>*{|fl|)  “  1  ”  ~  •  For  this  reason  we  introduce 

n 

the  following  two  function  spaces 


2  t  __  t 

F1  -  {«  e  W  '  (fl)  n  C2  I  (0)  -  0,  1  -  1/n  <  *  (  |0|  >  <  1} 


(25) 


12  * 

Fj  •  {*  e  w  '  (fl)  n  c2  i  *  (0)  •  o)  . 

An  approximate  version  of  the  constraint  ♦  (|fi|)  “  1 
itself.  This  is  demonstrated  by  the  following  proposition, 
constraint  t  (0)  “  0  is  superfluous. 


is  embedded  in  the  functional 
Note  that  as  in  section  1  the 
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Proposition  3.1.  J  ■  inf  J_(p) 

peF^  n 

is  independent  of  i,  i  ■  0,  1,  2. 

This  will  be  proved  in  conjunction  with  the  following  existence  result. 

Theorem  3.1.  Let  J  be  ss  in  Proposition  3.1.  Then  J_  is  minimized  on  P^,  that  is 

n 

there  exists  p.  e  F.  such  that  J  -  J _ ( 4<0 )  *  Also  if  J  is  minimized  by  some  p  € 

n  ^ 

F2(F2  3  F1 ) ,  that  is  J  «  j_(p)  then  we  can  redefine  p  on  a  set  of  measure  zero  so 

n 

that  p  satisfies  the  following: 

i)  p  is  continuous 

ii)  There  are  nested  convex  curves  Y^«  i  “  1,...,7r  -  1  such  that  Y^  “  {  p  “  i/n} 
for  i  *  2  and 

Y_  -  3(4)  -  1  -  1} 

n-1  n 


where  {4/  «  1  -  l/n}  is  a  convex  set.  Consequently 

n-1 

iii)  p  is  harmonic  in  8  \  {  U  y.  } . 

i-1 

Proof  of  Theorem  3.1.  We  first  introduce  a  simple  way  of  altering  a  function  p  e  P2 

with  J_(p)  <  «•  in  such  a  way  that  J  is  reduced.  To  carry  this  out  first  note  that 
n  _ 

>j-(p)  <  •  implies  that  | A„|  /  0  so  that  sup  p  >  1  -  1/n.  Now  recall  the  remarks  suds 

subsequent  to  the  definition  of  Cp.  For  each  i  •  1,...,F-  1  let  cj^n'  3  "  1»2  be 

the  convex  sets  associated  there  with  t  “  i/n  (if  sup  p  “  1  -  1/TT  then  we  only  have 

c!1*).  Let  p  be  the  continuous  function  which  is  i/n  on  r,  *  Jl11  c|^), 

n 

i  “  1,...,TC-  1,  and  is  harmonic  elsewhere.  Since  is  left  continuous  we  have 

V_(i/n)  -  V  (i/n)  and  therefore  j_(p)  <  J_(p).  By  [3]  we  see  that  p  e  Cj  so  p  e  F2 
as  well.  A  consequence  of  the  above  is  that  if  p  is  a  minimizer  then  p  -  p  a.e.  . 

This  establishes  the  second  part  of  Theorem  3.1. 

We  now  apply  the  above  methods  to  proving  existence.  It  is  easy  to  see  that  there 

exists  a  function  pQ  such  that  J_  (pQ)  <  »  so  choose  a  minimizing  sequence  (Pn) .  As 

—  n  _  ^  2 

above  replace  p  by  p  .  Choose  a  subsequence  such  that  p  +  p  weakly  in  W  '  (ft). 


As  in  the  proof  of  proposition  2.1  we  can  choose  a  subsequence  so  that  the  curves  yjn), 
associated  as  above  with  i|>n,  converge  uniformly  to  curves  y^  with 
i If  <  i/n  a.e.  outside  y^  and  tfi  >  1/n  a.e.  inside  y^.  Let  1)1  be  the  continuous 

function  which  is  i/n  on  y^  and  is  harmonic  elsewhere.  It  is  clear  from  the  uniform 
convergence  of  the  y*n*  that  1 ( 4>n )  J  ♦  |a^(i|*)|.  Also 

/  l??|2  <  J  |?*|2  <  lim  J  IV^I2 

o 

so  f  minimizes  J_.  Clearly  e  V  '  (ft) ,  ”  0  and  ^*(|ft|)  »  1  -  1/n  so 

n 

♦  e  f,. 


Proof  of  Proposition  3.1. 


Clearly 


J_  >  inf  J_ 
n  Fj  n 


since  F0  C  p.  C  f2.  From  theorem  3.1  we  know  that  there  exists  a  minimizer  of  J_  in 

n 

F^  and  that  if  i|i  e  F2  minimizes  J_  over  r2  then  <fi  e  F^ .  Therefore 

inf  J_  -  inf  J_  . 

F^  n  Fj  n 


TO  finish  we  will  show  that 


inf 


J__  <  inf  J_ 
n  F,  n 


by  altering  the  minimizer  +  whose  existence  was  proved  in  theorem  3.1. 

Let  C  be  the  interior  of  the  convex  set  {i|i  -  1-1/n}.  Recall  |c|  “  I An |  >  0. 
Since  the  Newtonian  capacity  of  a  point  is  zero  it  is  easy  to  construct  functions 


fn  e  CQ(C)  with  sup  fn  ”  1/rT  such  that 

/  |Vfn|2  -  o(n-’) 

and  |{f  «  1/n} |  «  0.  Now  consider  ^  +  fn*  claarly  6  po  and 


so  that 


J  <*  )  -  J  0|/)  +  o(  1/n) 
“■  n  — 
n  n 


Lim  J  (4>  )  -  J  <\|0 
—  n  — 
n+®  n  n 


and  we  are  done# 


-21- 


4.  Free  Boundary  Conditions  and  Applications 

Having  established  existence  for  the  approximate  problem  (23,24)  we  now  go  on  to 

establish  the  jump  conditions  across  the  free  boundaries  y ^  (see  theorem  3.1).  These 

will  be  applied  to  show  that  miniraizers  i>  of  J_  are  superharmonic  and  satisfy  L° 

n 


•  gradient  bounds. 

The  derivation  of  the  free  boundary  conditions  which  we  give  here  is  formal.  Our 
first  purpose  is  to  shed  light  on  their  consequences  so  we  have  left  the  justification  of 
these  formulas  to  another  paper  [9] . 

|  we  will  now  formally  apply  the  Hadamard  variational  principle.  The  result  will  be  a 

formula  for  how  the  normal  derivative  of  jumps  across  each  free  boundary  y^. 

Given  i  we  restrict  our  variations  to  *  ij  u  Ai+i-  Let  T£  be  a  one  parameter 
family  of  smooth  dif feomorphisms  of  onto  0^  which  coincide  with  the  identity  map  on 

some  fixed  neighborhood  of  3(1^.  Let  us  denote 


I 

Using  the  trial  function  <(>oT  we  compute  the  variation  of  J  . 

e  rT 

~  (J  (\|>oT  ))  -  /  |V*|2V  •  f  +  Vi|tT  •  pf  •  Vu 

dE  n  E  'e-0  fli 


-2 


(A, 


-  A 


i+1 


I  ~2>f  •  n 


where  jq  is  the  unit  normal  to  y^  directed  into  A^.  Hie  first  term  on  the  right  hand 
side  may  as  usual  be  turned  into  a  divergence  and  then  a  surface  integral. 


One  of  t..e  main  reasons  this  calculation  is  formal  is  that  it  must  be  justified  that 

variations  preserving  the  convex  function  class  W  n  Cj  provide  a  sufficiently  large 

collection  of  functions  f  •  n.  on  Y^  to  conclude  (26).  This  and  other  considerations 

such  as  smoothness  of  Y^  will  be  dealt  with  in  [9] . 

We  now  present  the  applications  of  (26). 

Theorem  4.1.  If  i(i  is  a  minimiser  of  J_  on  F  ^  then  i|i  is  superharmonic. 

n 

Theorem  4.2.  If  *  is  a  mlnimizer  of  J_  on  Ff  then 

n 

15,1  ‘  TvT  ’ 


Remark.  Recall  that  | ' |  <  Jf^)  for  solutions  of  the  variational  problem  (8,9)  and 

that  n  Vl A_|  is  the  difference  quotient  corresponding  to  ♦*,(|fl|).  From  this  and 
n 

other  considerations  we  expect  that  theorem  4.2  will  yield  an  estimate  independent  of  r\. 

As  a  precursor  to  theorems  4.1,  4.2  we  present  several  results  on  harmonic 
functions.  Let  A  be  an  open  region  bounded  by  two  nested  nonintersecting  convex 
curves.  Let  r ^  be  the  exterior  curve  and  the  interior  curve. 

Proposition  4.1.  If  h  is  subharmonic  in  A  with  h  -  on  i  •  1,2,c1  <  c2  then 

|Vh|  is  nondecreasing  along  curves  of  steepest  ascent. 
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Remark.  A  consequence  of  this  result  is  that  h  restricted  to  a  curve  of  steepest 
descent  (parameterized  with  respect  to  arclength)  is  a  convex  function.  This  provides  the 


intuition  that  the  volume  function  Vh<t)  =  |{h<t}|  should  be  convex.  This  is  made 
rigorous  by  the  following  result. 

Proposition  4.2.  If  h  is  as  in  proposition  4. 1  then 

V£(t)  <  0  on  (c1(c2)  • 

From  this  we  conclude  that 


ci  +  c?  c.  ♦  c, 

I  (c,  <  h  <  2  2  }  |  >  |  {  - -  <  h  <  c2)  |  . 

Proof  of  Proposition  4. 1 .  First  note  that  all  level  curves  of  h  are  convex  curves 

[3].  Given  p0  e  A  we  can  assume  without  loss  of  generality  that  pg  “  (0,0)  and  that 

(x0,y0)  «  (0,1)  is  the  outward  normal  to  the  level  curve  of  h  passing  through  pQ. 

Also  we  can  assume  that  y  «  gt(x)  gives  a  local  parametrization  for  level  sets 

{h  ”  t} .  That  is 


h(x,gt(x) )  -  t 

for  x , t  small.  Note  that  gj(0)  -  0,  gg(0)  <  0.  From  the  chain  rule  we  get 

hxx  +  hxy*t  +  hyy<9£>2  +  Vt  *  0  * 

Using  Ah  >  0  and  g^(0)  »  0  we  have 


Also 


and 


“hyy  +  hygg  >0  at  p0  . 


Vh 

Tvhl 


•  7(|7h  r> 


2  I  Vh ! -1 (h2h  + 

'  X  XX 


2hxhyhXy  + 


h^yy) 


hx(p0)  =  0*  ~hy  -  |7h| 

SO 

fSr  *  V<,7h|2>  ’  "2Vyy  *  -2hygo 

at  Pg.  Since  level  curves  are  convex  and  (0,1)  is  an  outward  normal  we  see  that 

gg ( 0 )  <  0.  Therefore  |Vh|  decreases  as  we  go  outward  along  curves  of  steepest  descent. 
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Proof  of  Theorem  4.1.  Assume  p  is  a  minimizer.  To  prove  that  t|i  is  auperharmonic  we 
will  show  that 

/  7p  •  7  n  >  0 

a 


for  all  n  e  C*(8).  For  simplicity  we  prove  this  only  for  n  e  C^a^),  1  -  1, 

8^  -  u  Ai+1.  We  use  the  divergence  theorem  to  evaluate  the  above  integral.  Since  the 
free  boundary  condition  is  only  weakly  attained  the  following  calculation  is  rigorous  only 
with  the  addition  of  a  limiting  argument.  The  technical  details  of  such  an  argument  are 


left  to  [9] . 


Formally  the  divergence  theorem  says  that 


From  the  free  boundary  condition  we  see  that  the  difference  of  th  normal  derivatives  is 
positive  if  | A± |  >  |Ai+1|  and  negative  if  ) f  <  |a1+1|.  In  the  first  case  would 
be  superharmonic  in  8^  and  in  the  second  subharmonic.  Assume  | A^ (p )  |  <  |  A^+1  (i|i  )  |  so 
that  \|i  is  subharmonic  in  Bi .  Let  h  be  the  harmonic  function  on  8i  which  agrees 
with  ip  on  y  and  Yi+1-  Clearly  ♦  <  h  in  8^  From  proposition  4.2  we  have 
|  A± ( h )  I  >  |A1+1(h)|  but  since  <  h  the  level  set  {\|»  »  i/n}  must  be  surrounded  by 

{h  =  i/n}  in  which  case  |Ai<i(i)|  >  |Ai+1<$>|.  But  this  contradicts  our  assumption  so  we 
conclude  that  |Ai(p)|  >  and  80  ♦  is  superharmonic. 

Proof  of  Theorem  4.2.  To  establish  the  gradient  bound  we  use  proposition  4.1  and  the  free 
boundary  condition.  Given  x  e  for  some  i,  consider  the  curve  of  steepest  ascent 
which  leads  from  x  to  3An-  By  proposition  4.1,  1 7ip |  is  nondecreasing  along  this 

curve  except  across  the  free  boundaries  where  the  decrease  is  given  by  the  free  boundary 
condition.  Thus  we  have 


as  required 


|7t(x)|2  <  ^  n  "1{|^+1 
k«i 


-2 


lAk|-2) 


<  n  -2/|^  r2 


Remark .  Ae  results  of  this  paper  except  theorem  2.2  have  been  generalized  to  n- 
dimensions.  This  will  be  dealt  with  to  some  extent  in  future  papers.  In  addition  soeie  of 
the  results  generalize  to  the  variational  problems  associated  with  functionals  of  the  type 


/  |V*|2  +  (Y  -  I)"1  /  U(*)(***)Y,  Y  >  1 

a  o 


where  u  is  smooth  and  satisfies  u  >  X  >  0.  This  more  general  class  is  used  in  [8]  to 
model  adiabatic  compression  of  a  plasma. 
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